AD=A102 665 NAVAL POSTGRADUATE SCHOOL MONTEREY CA F/6 11/6 
AN ELASTIC@PLASTIC FINITE ELEMENT ANALYSIS OF NOTCHED ALUMINUM =-ETC(U) 
MAR 81 M J KAISER : 

INCLASSIF IED 


LS 
COE 


EEEEEEEEER 

Pe 
BREE EEEEEeee 

EREEEEEREREEES 


LEVEL! fy 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


Te) 
Ne) 
Ne) 
QL 
—) 
= 
= 
[am 
<< 


THESIS 


AN ELASTIC-PLASTIC FINITE ELEMENT ANALYSIS 
OF NOTCHED ALUMINUM PANELS 


by 


Michael John Kaiser 


March 1981 


Thesis Advisor: G. H. Lindsey 


Approved for public release; distribution unlimited 


UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE (When Dave Entoreg) 


READ INSTRUCTIONS 
DRE COMPLETING FORM 


2. -Abd FUORI 0} - 
Michael John| Kaiser 


9. PERFORMING ORGANIZATION NAME AND ADORESS 0. PROG 4 ELEMENT. PROJECT, TASK 


REA ORK UNIT NUMBERS 
Naval Postgraduate School 
Monterey, California 93940 


11. CONTROLLING OFFICE NAME ANO ADORESS 


Naval Postgraduate School 
Monterey, California 93940 


Init sista. cio i chip icin i tel ie ik ie 


Approved for public release; distribution unlimited 
17. OISTRIGUTION STATEMENT (of the ebetract entered in Block 20, !f different trem Repert) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree side if nececcary and identify by bieck number) 
Finite element analysis, stress, stress concentration factor, 
notched panels, residual stress, plastic zone, ADINA,extrapolatio 
elastic-plastic stress. 


20. AQSTRACT (Cantinue on reverse cide if necesoary and identify by bieck mumber) 

Finite element, elastic and plastic analyses of various alumi- 
num panels, containing holes and notches, were conducted for com- 
i parison with photoelastic experimental results. A FORTRAN IV 

| program, ADINA (Automatic Dynamic Incremental Nonlinear Analysis), 
was used for both linear and nonlinear analyses. Mesh refinements 

were used for each panel and the monotonically convergent results 
were extrapolated using Richardson's method. Stresses were .. 


{ DOD an 72 1473 sence Ne ' MOV 6815 OBSOLETE | rave UNCLASSIFIED 
, Lhah vier P eeDRTT CLR ETON OF FE FLOR Te Bae ee 


25 4S OC 


UNCLASSIFIED 


ae 


Item 20 (contd). 


ocally smoothed from the Gauss integration points to the nodal 
points. Eight noded, isoparametric elements were used throughout. 
Modification to an ADINA preprocessor program, also coded in 
FORTRAN IV, was made for use with a VERSATEC plotter. 


Comparisons were made to the elastic, analytic series solution 
by Howland for a circular hole in a finite strip. The finite ele- 
ment results varied by less than one percent from Howland's solu- 
tion. Handbook values for the elastic stress concentration factors 
of the geometries investigated differ from finite element results 
s by less than one percent in all cases. The photoelastic works of 
Frocht were also compared where applicable. g”yStresses in the plasti 
range obtained from slip-line theory for a rliigid-perfectly-plastic 
material show excellent correlation to a finite element analysis of 
such a material. Comparisons to elastic ang plastic experimental 
data were made for the panels analyzed and ghow good correlation 
to finite element results. 


sion For 
CRA&L 

perc TA® = 

Unannenoced ro 


4 orn. —---~ 
Just ypicet tom 


Acces 
“NIIS 


_—_—" 
_— 


asl srmaget iel 


DD tay 1473 j 
S/ N 104-014-6601 SECUMETY CLABNMEATION OF THIS PAEEEDER Bere Enreves) | 


Approved for public release; distribution unlimited 


An Elastic~Plastic Finite Element Analysis 
of Notched Aluminum Panels 


by 


Michael John Kaiser 
Lieutenant Commander, United States Navy 
B.S., St. Cloud State University, 1969 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1981 


Author: 


Approved by: 
Thesis Advisor 


tU. F. Fee 


ke_a_Few—n nae ee — 
Chairman, Derartment of Aeronautics 


Dean of Science and Engineering 


ABSTRACT 


Finite element, elastic and plastic analyses of various 


aluminum panels, containing holes and notches, were conducted 
for comparison with photoelastic experimental results. A 
FORTRAN IV program, ADINA (Automatic Dynamic Incremental 
Nonlinear Analysis), was used for both linear and nonlinear 
analyses. Mesh refinements were used for each panel and the 
monotonically convergent results were extrapolated using 
Richardson's method. Stresses were locally smoothed from 
the Gauss integration points to the nodal points. Eight 
noded, isoparametric elements were used throughout. Modi- 
fication to an ADINA preprocessor program, also coded in 
FORTRAN IV, was made for use with a VERSATEC plotter. 

Comparisons were made to the elastic, analytic series 
solution by Howland for a circular hole in a finite strip. 
The finite element results varied by less than one percent 
from Howland's solution. Handbook values for the elastic 
stress concentration factors of the geometries investigated 
differ from finite element results by less than one percent 
in all cases. The photoelastic works of Frocht were also 
compared where applicable. Stresses in the plastic range 
obtained from slip-line theory for a rigid-perfectly-plastic 
material show excellent correlation to a finite element 
analysis of such a material. Comparisons to elastic and 
plastic experimental data were made for the panels 


analyzed and show good correlation to finite element results. 
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I. INTRODUCTION 


Development of on-board fatigue monitoring systems 
for Naval aircraft have made it possible to record exten- 
sive structural loading data in flight. The strain gages 
used in such a system must be located away from stress 
concentration areas to prevent their fatigue; however, 
these areas are of the greatest interest in analyzing 
and predicting fatigue life of the structure. Understand- 
ing the relationship between nominal, far-field stresses 
and local stresses in critical areas thus becomes vitally 
important. Recent experimental investigations into the 
effect of uniform, far-field loads on stress concentration 
areas have been made at the Naval Postgraduate School 
(NPS) using photoelastic techniques [Refs. 1, 2 and 3]. 
These experiments involved loading 7075-T6 aluminum into 
both the elastic and plastic regions, as well as measure- 
ments of residual stresses resulting from plastic yielding. 

Finite element analyses (FEA) of the aluminum panels 
used in the experiments of Stenstrom [Ref. 1] were con- 
ducted. The panels used in the experiments of Engle 
[Ref. 2] and Stuart [Ref. 3] have similar geometry. The 
finite element programs available at NPS were surveyed 
and ADINA [Ref. 4] was chosen for its proven ability to 


produce the nonlinear analyses required for plastic 
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yielding of aluminum. To provide increased accuracy, 
each panel was modeled using two meshes. The results 
obtained for the coarse and fine meshes were extrapolated 
to a final result using the Richardson extrapolation 
technique [Refs. 5 and 6]. 

Along with ADINA, a preprocessor program, PSAP1 [Ref. 7], 
was used to verify mesh connectivity prior to analysis by 
ADINA. PSAP1 provides a graphical output of the finite 
element mesh and was coded for the CALCOMP plotter installed 
at NPS prior to 1978. For this thesis PSAPl was adapted 
for use with the VERSATEC plotter now installed at NPS. 

The stress-strain material properties of the 7075-T6 
aluminum actually used to make the panels had to be estab- 
lished to provide an accurate material model for use with 
ADINA. Material testing was conducted to establish the 
Young's Modulus (E), Poisson's Ratio (v), yield stress (a), 


strain hardening modulus (E,) and the Ramberg-Osgood co- 


fawn Kacibet sence et tie poten iy hh Dhaene tara Ralltierenents aliens angmntemanin ot 


efficients 8 and n. 

Comparisons were made to other works, in addition to 
the experiments conducted at NPS. The initial analysis 
involved a comparison of FEA to the results of Howland 
[Ref. 8], for a circular hole in a finite strip, to vali- 
date the methods used. A comparison of FEA to plane 
stress, Slip-line theory, for rigid-perfectly-plastic 
material was also included as a validation for the plastic 


analyses. 
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II. MATERIAL PROPERTY TESTING OF 7075-T6 ALUMINUM 


The elastic and plastic material properties of the ; 
aluminum panels were established by tensile tests of uni- 
axial specimens made from the same mill run. The speci- 
mens were manufactured and tested according to current 
ASTM standards [Ref. 9]. MICRO-MEASUREMENTS, EA~-13-125AD- 
120, precision strain gages with a temperature compensated 
bridge circuit were used on all specimens. Transverse 
gage sensitivity errors were corrected according to the 
manufacturer's recommendations [Ref. 10]. Critical cross- 


section measurements were made with a micrometer. 


A. TESTS FOR AXTAL LOADING 

Initial tests of the two-gaged specimen, Fig. 1, in 
the MTS testing machine indicated a significant bending 
moment was being produced by the 30,000 1b GRIFF grips. 
To investigate this problem further, tests were conducted 
on both the MTS and RIEHLE test machines with a five- 
gaged specimen shown in Fig. 2. The results of these 
axial loading tests, shown in Table I, verified that the 
GRIFF grips on the MTS test machine do not give axial 
loading. An inspectionof the gripped region on the speci- 
men indicated that the jaws of the grip were not applying 


a uniformly distributed force and thereby induced a bending 


3 


moment by off-axis loading as shown in Fig. 2. The grips 
on the RIEHLE test machine gripped evenly and a uniform 


strain distribution resulted as seen in Table I. 


B. CHARACTERISTICS OF 7075-T6 ALUMINUM PANELS 
The following characteristic properties were determined 
from the four specimens tested. 
1. Young's Modulus (E) 
Tests were conducted using the specimen shown in 
Fig. 3 on the MTS test machine with 10,000 1b INSTRON grips, 


which gripped the specimen evenly. The results of testing 


one 


three specimens are shown in Tables II to IV. Linear 
regression in the elastic range of all the test results 


6 psi, with a 


determined a Young's Modulus of 10.12 x 10 
correlation coefficient of 0.9996. 
2. Poisson's Ratio (v) 


Tests were conducted using the specimen shown in 


} Fig. 1 on the RIEHLE test machine with 10,000 1b RIEHLE 


Re DE ee a a a a ee ee eee eee, 


: grips. The results are tabulated in Table V. Linear 


regression of these results in the elastic region deter- 


mined Poisson's Ratio to be 0.3256 with a correlation if 


coefficient of 0.99996. 


3. Yield Stress and Strain Hardening Modulus 


These values, required for ADINA's bi-linear 
material model, were determined graphically using the 


data from Tables III and IV. Plastic region data in 
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Table II is not reliable because of excessive creep en- 


countered during that test. 


0.2% offset yield stress, o, 76,000 psi 


strain hardening modulus, E. = 566,000 psi 
The graphical fit of these values to the test data can be 
seen in Fig. 4. 
4. Ramberg-Osgood Coefficients 
The Ramberg-Osgood equation for elastic-plastic 


Seem noni tinal mnedeaibermenicediiaipe, 


stress-strain characterization is given by: : 


_ oO o,n 
e=5 + 8(5) (1) 


where: 
€ = strain 
ag = stress 
E = Young's modulus. 
The 8 and n coefficients were determined graphically from 
} the data of Table IV, by the method given by Rivello [Ref. 11]. 
- The data in Table IV gave the following values which are 
the best fit to the combined test data 


43 


8 = 1.479 x lO . 


a n = 21.58 


The graphical fit of these values, in Eq. (1), with the 


test data is also shown in Fig. 4. 


| 
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MODIFICATION TO GRAPHICAL PREPROCESSOR 


The use of a graphical preprocessor program, such as 


PSAP1, is vital in detecting mesh errors that may other- 


wise go unnoticed. Establishing the proper node locations 


and element geometry prior to analysis for a complex code 


such as ADINA is of utmost importance. 


A. PSAP1 MODIFICATIONS 


The program PSAP1 was originally coded in FORTRAN IV 


for use on the NPS IBM 360/370 installation with the 


CALCOMP Model 765 drum plotter. The CALCOMP system o: 


installed at NPS used the +Y axis as the unlimited plotting 


axis, see Fig. 5. The entire plotting logic in PSAP1 uses 


this orientation of axes to allow multiple plots in a con- 


tinuous strip. With the VERSATEC Model 8222A electrostatic 


plotter now installed at NPS, the +X axis becomes the un- 


limited plotting axis, shown in Fig. 5. To avoid an ex- 


tensive recoding of PSAP1 for use with the VERSATEC plotter, 


a Simple coordinate transformation of the plot was made in 


a limited number of short subroutines. First, all installa- 


tion dependent plotting calls used in PSAP1 were identified. 


These involved seven plotter functions for which new sub- 


routines were coded. 


4 
: 


Function New Subroutine 


SS i li a a ll et SS 


Initialize Plotter CALCMP 
Move Plotter Pen CALPLT 
Letter on Plot NOTATE 
Number on Plot CALNUM 
Determine Current Pen Location CALWH 
Stop Plotter PSTOP 


| 
| 
Draw a Line CALINE 
The subroutines listed above merely rotates the plot 

to coincide with the VERSATEC axis orientation and retain 

all features originally in PSAP1. Since all plotter | 
hardware code is now isolated in these seven subroutines, 

future adaptations to other plotting systems is simplified. 


To provide documentation of this update to PSAPl, a complete 


listing of the new program is provided in Appendix D. 


B. USE OF PSAPI 

Previous use of PSAPl1 on the IBM 360 system necessitated 
use of a load module since PSAP1 took over one minute of 
CPU time to compile. With the new IBM 3033 system 
compilation requires eight seconds; however, use of a load 
module or disk stored source code is still recommended 
since PSAP1 contains roughly 2,500 lines of code. Appendix 
A contains sample JCL to use PSAP1 on the IBM 3033 MVS 


system. With minimal effort PSAP1 could also be set up 
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for use on the IBM VM/370 system. The user's manual for 


PSAP1l is in Ref. 7. In addition to a mesh plot, PSAP1 


provides a listing of the node coordinates, element 


connectivity and several key input values used in execu- 
tion of ADINA. This information provides a useful check 


of the input data. 
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IV. FINITE ELEMENT ANALYSIS (FEA) 


| A. DESCRIPTION OF MESHES USED 


1. Length to Width Ratios and Boundary Conditions 


Initial finite element models of specimens had 


length to width ratios near one, as used by Garske [Ref. 12] 
for his FEA, but they did not provide the desired uniform 
distribution at the loading boundary. Specimen length to 
width ratios of 3-5 were used by Armen, Pifko and Levine 
{Ref. 13] in their FEA and by Stenstrom [Ref. 1] in his 
photoelastic experiments. The criteria established to 
determine uniform boundary stress distribution was uni- 
formity in nodal displacements along the loaded edge as 
discussed by Segerlind [Ref. 14]. In the models used for 
FEA in this thesis, nodal displacements were uniform to 

‘ within 0.1%, and the resulting stress distribution was 
uniform axially to within 0.1% at the panel ends. In all 
cases two-dimensional, eight noded, isoparametric elements 


were used. These higher order elements cannot be loaded 


in an "intuitive" manner as discussed by Zienkiewicz 
{Ref. 15, p. 223]. Figure 6 shows the nodal loading 


required to obtain a uniform surface load. 
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2. Element Mesnes 
Two meshes were developed for each panel analyzed. 
A reasonable effort was made to keep element corner angles 
as close to 90° as possible to reduce the effect of element 
distortion discussed by Hopkins and Gifford [Ref. 16]. 


All meshes modeled a quarter of the actual panel by using 


TS dedi oe 


the two axes of symmetry as is common practice in FEA. 

The step from course to fine element meshes was made so 
that each element in the course mesh was subdivided into 
four smaller elements of the same type. Such a mesh sub- 
division can be expected to give monotonic convergence of 
results, Cook [Ref. 17], and allow extrapolation to results 
of an infinitely fine mesh. Figures 7 through 13 illus- 
trate the element meshes used in this analysis as plotted 


by PSAP1. 


B. COMPUTATIONAL PROCEDURES 

1. Using ADINA 

Once the mesh has been developed, input data is 

prepared in accordance with the ADINA user's manual [Ref. 4]. 
This same set of data is then used as input for PSAP1 to ~ 4 
check for errors and provide a graphical display of the 
element mesh. After preprocessing by PSAPI1, the data is 
entered into ADINA for analysis. In the case of linear 
analysis, two types of stress output may be specified, 


nodal point or Gauss integration point. Nodal point output 
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can be computed for up to eight node point stresses for 

each element. Since 2x2 Gauss integration was used, four 
Gauss point stresses were computed for each element. The 

2x2 Gauss integration is recognized as the most efficient 
integration order for this type of analysis [Ref. 15, p. 284]. 
The linear analysis used an isotropic linear elastic material 
model (MODEL "1" in Ref. 4) which required input of E and 

Vv material properties. The nonlinear analysis allows only 
Gauss point stress outputs and uses a bilinear elastic- 
plastic material model, with von Mises yield condition and 
isotropic strain hardening (MODEL "8" in Ref. 4). 

For static analyses ADINA uses a time function 
method to apply loads in steps. Linear analysis loading 
was accomplished in a single step to a nominal value of 
3,000 lbs load. Nonlinear analysis loads were applied in 
ten steps to a maximum value, matching the experimental 
loads, and then unloaded to zero in ten steps to obtain 
residual stresses. The stress output from ADINA is a 
listing of nodal or Gauss point stresses for each element. 
Since the only area of interest in this analysis was the 
distribution of stresses along the reduced cross-section, 
no large post-processing program was developed or used. 
All final computations using ADINA output data were accom- 
plished on a HEWLETT-PACKARD 9830A calculator, using short 


programs coded in BASIC. If more extensive stress 
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distribution information were desired, some form of auto- 
mated post~processing would be necessary to reduce the 
computational workload. At a minimum, nodal stress outputs 
by ADINA must be averaged to obtain unique values of stress 
at nodes shared by more than one element. 
2. Richardson Extrapolation 

The use of course and fine meshes allows extrapola- 

tion to an infinitely fine mesh as discussed earlier. 


Richardson extrapolation [Ref. 5] was used in this analysis 


where: 
a (hi)™ = oi(h_)™ 4 
g = —< 5 aii (2) i 
extrap h™.y,™ ; 
B ¢ , 
iq 
where ‘ 
if 
= 1 iq 
) Coxtrap extrapolated solution 4 
oa = solution obtained with h=h, : 
on = solution obtained with h=h, ft 
f 
: he = linear dimension of course element : 
ha = linear dimension of fine element d 
\ 
m = 2 (for this analysis) * 
The exponent m is determined by the order of the discre- : 
tization error O(h"). Since h represents the length of an 1 


element the element area is represented h?, In a two 


dimensional problem such as this O(h™) is of the order of 


{ h?, the area of an element. In the mesh refinement scheme 
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h 
used h, = sh or = = 5: Equation (2) can be rewritten 
- c 
of Ned =o (Pe)? 
c ne F Ae 
a = (3) 
extrap hp ty 
ae = (-=) 
c c 
thus Ll 
Oo, - FO 
ees 4 “c 
Sextrap = E} 3) 
4 


Equation (4) then becomes the relation to obtain extrapolated 
stresses from coarse and fine mesh results in a two dimen- 
sional analysis. Better extrapolations can be obtained by 
using three or more refined meshes, but, the computational 
effort increases significantly. 

3. Optimal Stress Locations and Local Smoothing | 

It is generally accepted that the most accurate 

sampling points for stresses are the Gauss integration i 
points within the element [Ref. 15, p. 281, and Ref 18]. 
In this analysis, the nodal points are of the greatest i 
interest; thus a technique of local smoothing must be ; 
applied to the integration point stresses to obtain nodal ~ oH 
stresses as reported by Hinton and Campbell [Ref. 19]. 
The formulation of this local smoothing technique for ADINA 
elements is developed in Appendix B. The nodal values 4 
obtained must then be averaged if shared by two or more 


{ elements. ; 
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4. Computational Times 


Because of the extremely large size of ADINA (about 
17,000 lines of code in the NPS version) and the out of core 
solver, it does not adapt well to time sharing systems. 
Using the IBM 360 system at NPS, ADINA required 31 user 
defined overlays to create a manageable load module in 
about 30 minutes of CPU time. With the new IBM 3033 MVS 
system at NPS, ADINA is compiled without overlays in about 
one minute. When using a load module, the program execu- 
tion took less than 2 minutes CPU time on the IBM 3033. 

In addition to ADINA, the preprocessor (PSAP1) 
and post-processing techniques involve considerable time 
and effort. Figure 14 is a flow chart of the computational 
procedure used in this analysis. An example of the JCL 
to use ADINA in load module form on the IBM 3033 MVS 


system with use of the mass storage facilities is shown 


in Appendix C. 


re 


V. RESULTS OF ANALYSIS 


A. CIRCULAR HOLES IN LINEAR MATERIAL 

The FEA results for a circular hole in a finite width 
strip were used to validate the elastic computational pro- 
cedure discussed earlier. The results of Howland [Ref. 8] 
were compared to both the Gauss point smoothed results and 
the nodal output results in Fig. 15. The stress concentra- 
tion factor o/s, is referenced to the far-field stress. 
The smoothed results give the best match to the results 
of Howland at the edge of the hole, and the only signifi- 
cant variation between the two FEA methods occurs within 
the first 0.25 inches from the edge. In order to obtain 
the 0.25 inch stress value for the coarse mesh, in the 
Gauss point smoothed result, a midside node value had to 
be obtained by the averaging method discussed in Appendix 
B. The linear distribution of smoothed stresses alcng the 
sides of the element, [Ref. 19], appears to produce a 
less accurate result in this area of extreme stress 
gradient, when compared to ADINA's nodal output result. 
This tendency was noted in all cases; however, the peak 
stress values from smoothed results consistently gave 
better correlation with other investigators [Ref. 20]. 

A circular hole with A=0.25 was also analyzed and 


compared to the experimental results of Stenstrom [Ref. 1] 
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along with an interpolation of Howland's results. The Bg 
experimental data correlates well with the FEA results; 
however, the a. experimental data shows significant varia- 
tion between 0.125 and 0.375 inches from the edge of the 


hole, as seen in Fig. 16. 


B. OPPOSITE U NOTCHES IN LINEAR MATERIAL 
The results for linear analysis are presented in non- 
dimensional stress concentration form; however, the 
normalizing stress changes. For U notches Kp is the 
stress concentration factor referenced to the theoretical, 
nominal stress (o.) in the reduced cross-section where 
oS Load/Area of Reduced Cross-Section. 
1. Shallow Notch Panel 
The FEA results plotted with the experimental 
data of Stenstrom are shown in Fig. 17. Once again both 
FEA results are shown and the variation for the two methods 
occur within 0.25 inches from the notch edge; however, 
there was less variation than was seen in the circular 
hole analyses. 
For this panel the experimental data appear to be 


uniformly below the FEA results for o The OY data shows 


g 
significant variation at the 0.625 inch point but follows 
the proper trend within 0.5 inches from the notch edge. 
According to data collected by Peterson (Ref. 20], this 


notch geometry should yield a maximum Kn = 2.74. The 
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Gauss point smoothed results matched this value exactly. 


Stuart [Ref. 3] reported a K, = 2.69 for this same notch 


geometry, with a standard deviation of 0.187 in the 14 
samples he measured by photoelastic methods. Early photo- 
elastic work by Frocht [Ref. 21] determined Ky, = 2.7 for 
this notch geometry. The FEA results appear to be in good 


agreement with other investigators, for this notch geometry. 


2. Deep Notch Panel 
Results of the FEA for a deep notch were plotted 


with Stenstrom's experimental data in Fig. 18. The results 
of the two FEA methods again diverge within 0.5 inches from 
the notch edge. FEA stress values at the 0.25 inch point 
have spread farther apart in this case since the stress 
gradient is very severe at that point. The experimental 


data correlates well for both o, and o.) however, the maxi- 


9 


mum experimental Tg is considerably lower with a Kn = 3.83. 


Results reported by Stuart for this notch geometry was a 


Kn = 4.05 with a standard deviation of 0.219 for 14 speci- 


mens measured by photoelastic methods. Frocht [Ref. 22] 


reported a photoelastic K,, = 3.9 for this notch geometry, 


T 
but concluded that the result was 5-10% low, giving a 


corrected range of K,, from 4.1 to 4.3. The Gauss smoothed 


T 


FEA result gave a K,, = 4.24 which compares well to an 


T 
empirical relation given by Peterson [Ref. 20] for 


r/d < 0.25. 


2t 


toc aoe 2t, 2 2t,3 | |. 
Kya (1 - $=) (0.78 + 2.243vt7x) 0.993 + 0.18 = - 1.06() + 1.7115) (5) 


where 
t = notch depth (3.9375) 
r = notch radius (0.625) 
d = minimum width (15.625) 
D = maximum width (23.5) 


Inserting the values above into Eq. (5) gives a K, = 4.26, 


T 
which is 1/2% above the FEA result. Other FEA results 
reported by Armen, Pifko and Levin [Ref. 13] and Griffis 
[Ref. 23], using linear strain triangle (LST) elements, 
produced Kn values within 5% of those produced by use of | 
Eq. (5). The rectangular elements used in this analysis 

are known to give better results than LST elements as 

noted by Clough [Ref. 24]. It is clear that the FEA 


results obtained for this notch are in good agreement 


with other works. 


C. OPPOSITE U NOTCHES IN NONLINEAR MATERIAL 

The analysis for loading into the plastic region of the 
7075-T6 aluminum was made using the bilinear material 
model discussed earlier. The loads used were selected to 
match those used in the experiments of Stenstrom; thus 
allowing direct comparison. The strains obtained in those 
experiments were used to solve for stresses by use of the 


Prandtl-Reuss plastic flow equations. 


ee ee 
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1. Shallow Notch Panel 
The results of FEA for the three load cases, 60,000, 
65,000 and 70,000 lbs are presented along with the experi- 
mental results in Figs. 19 through 21. The oi results 
compare well although no trend for peak 4 stress away from 
the notch edge is shown in the experimental data. In all 
cases the FEA determined the peak J, stress to occur near 


the yield boundary, and the gradient of the o, stress to 


8 
fall off dramatically in the plastic zone. This character- 
istic behavior of the J, stress was reported by other 
investigators (Refs. 13 and 23] using FEA on 2024-T3 aluminum. 
Plane elastic-plastic stress distributions reported by 
Frocht [Ref. 25] show similar trends. The experimental 

data also shows a marked change in the gradient of 5, stress 
within the plastic region. The growth of this plastic 
region is approximated using the FEA results for this notch 
in Figs. 22 through 24. 

Experimental data for the J, stress distribution 
matches the FEA results closely except at the notch edge 
where the measured C, does not go to zero as it should. 

The characteristic peak value of d. near the plastic 
boundary as seen in the FEA results is also shown by the 
test data. 

The FEA residual stress computed upon unloading 


from the three load cases are shown in Figs. 25 through 30. 


The characteristic distributions of the Tg residual stress 


agrees with those reported by others [Refs. 25, 26 and 27]. 
The experimental residual stress distributions reported by 
Stenstrom show similar trends but significant variations 
when compared to the FEA results. 
2. Deep Notch Panel 

Three load cases were computed to plastic loading 
levels; however, only limited experimental results are 
available for this notch as seen in Fig. 31. The 30,000 lb 
load is just at the onset of yield in the notch root area. 
Limited residual experimental data [Ref. 3] was available 
for comparison in Figs. 32 and 33 which are plots of the 
residual Ty and OL stress distributions as a result of the 
three loading cases. The 30,000 lb load has caused yielding 
in a small region at the root of the notch as seen in Fig. 34. 
Figures 35 and 36 illustrate that the plastic zone does not 
grow to the extent it did in the shallow notch. The stress 
gradients are very severe in the deep notch and as a result 
the 0.25 inch sampling points used in the FEA may be useful 
in only showing gross trends close to the notch edge. The 
trends appear to be much the same as in the shallow notch, 


with the peak g. and o 


3 r stresses occurring near the yield 


boundary. The yield boundaries shown in Fig. 31 are 


approximations based on qualitative analysis of the finite 


element results. The Ty experimental data for the 30,000 


lb load case correlates well with the peak stress, again 
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appearing low as it did in the linear analysis. The residual 


Oo, Stress distribution in Fig. 32 follow much the same trends 
as seen in the shallow notch but with much higher gradients 
within the first 0.5 inches from the notch. The FEA results 
for the residual ona stresses shown in Fig. 33 indicate a 
limitation in the element size used since the Oe value of 
the notch edge does not return to zero as it should. Because 
of this problem the data may be questionable for showing 
proper trends in the first 0.5 inch from the notch edge. 
3. Rigid-Perfectly Plastic Panel 

A stress distribution for the theoretical material 
used in slip-line theory was desired. By using ADINA's 
bilinear material model such a material could be approxi- 
mated reasonably weli. For this analysis a Young's Modulus 


(E) of 1076 


psi was used to model perfect rigidity. The 
strain hardening modulus (E,) was set to zero to model 
perfect plasticity. Poisson's ratio (vv) was 0.4999999, 

as close to 0.5 as the computer would allow. The E = 1076 
also represents a computer limit in approaching an infinitely 
large E. Figure 37 illustrates the results obtained and 
compares them to a slip-line solution. The results are 
normalized to the arbitrary 73,000 psi yield stress used 

in the analysis. The Ty values obtained agree exactly 

with slip-line theory. Conversely the oy results do not 


reflect the same values as slip-line theory, but do show 


a similar trend. The growth of the plastic zone obtained 


is shown in Figs. 38 through 40. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The results obtained from the FEA have proven useful 
in determining the validity of experimental data gathered 
by photoelastic techniques. The FEA results varied by 
less than 1% when compared to published analytical results 
and handbook values. Experimental data from Stenstrom's 
photoelastic work correlated well with the FEA results. 
The primary exception was the residual stress experimental 
values, which varied significantly from the FEA results. 
The variation may be in transforming the residual strains 
measured photoelastically into stresses for comparison 
with the FEA results, since ADINA only provides a stress 
output. Limitations of ADINA's bilinear material model, 
initially considered severe, do not appear to have hampered 
this investigation. A possible exception is the residual 
analyses where the transition region from elastic to 
plastic strains becomes especially important. The Gauss 
point smoothed results gave the best correlations at the 
edge singularities in all cases; however, due to the 
limitations noted at 0.25 inches from the edge, the 
results at that specific point may not be as accurate for 
this method. The nodal output results gave consistently 


higher stress values at the edge singularities. Because 
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of the severe stress gradients near the deep notch analyzed 
the use of a finer element mesh near the notch would probably 
produce better results. 

The effort involved in developing two meshes such as 
those used in this thesis is considerable. An automatic 
mesh generation capability would reduce the workload and 
allow experimentation with several types of element meshes. 

ADINA proved to be a useful and powerful program, as 
expected, but something simpler and less awkward to use may 
be all that is required for two dimensional analysis. Such 
a system is already in use at NPS but does not offer non-~ 
linear capabilities. If use of ADINA is to be continued in 
this type of investigation, a post-processing program should 
be adapted. There are programs available to post-process 
ADINA data at NPS [Refs. 28 and 29] but they would require 
modifications to work with two dimensional analyses and 
the VERSATEC plotter. 

Standardized material property testing would ease the 
inevitable task of obtaining basic material properties for 
use in analysis or experiments. Some form of automatic 
data collection with use of the MTS testing machine would 


allow testing of a larger sample population and provide 


statistically more accurate information. 
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TABLE’ I. MTS AND REIHLE 5 GAGE TEST RESULTS 
All Strains are 107° in/in 


MTS Test Machine 


! LOAD STRAINS 
lbs ey e5 e3 Ey Es 


REIHLE Test Machine 


Nee ete ae lle oo 


LOAD STRAINS 
lbs 


1564 1540 


2000 
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TABLE II. MTS SPECIMEN A TEST RESULTS 


Cross-section = 0.03975 in 


Load, lbs. Strain 10— 
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TABLE IV. MTS SPECIMEN C TEST RESULTS 
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TABLE VI. 


TABLE VII. 


A = 0.2 HOWLANL DATA 


DISTANCE 
FROM HOLE 


X = 0.25 HOWLAND 


DISTANCE 
FROM HOLE 
in. 


TABLE VIII. A = 0.2 FEA RESULTS - NODAL 


o, = 5000 psi 
DISTANCE 
FROM HOLE, in. 


} = 0.2 FEA RESULTS ~- GAUSS OUTPUT 


5. = 5000 psi 


DISTANCE 
FROM HOLE, 


TABLE X. d = 0.25 FEA RESULTS - NODAL OUTPUT 
DISTANCE ABOR CPSs 
FROM HOLE, in. 


TABLE XI. 0.25 FEA RESULTS - GAUSS OUTPUT 


o., = 4500 psi 


DISTANCE 
FROM HOLE, in. 0, ,psi 


TABLE XII. SHALLOW NOTCH FEA LINEAR RESULTS - NODAL 


a 4225 psi 


DISTANCE 
FROM NOTCH, in. o, ,psi Oo. ,psi 


TABLE XIII. SHALLOW NOTCH FEA LINEAR RESULTS - GAUSS 


oo 4225 psi 


DISTANCE 
FROM NOTCH, in. Ty *psL ee *PS1L 
fia 
4.25 Srar.d Lif?agr 
4,568 | Pasoes 1284.3 be 
a.7s Sonn, 9 ris.s 
1.468 Seal. 2 Prana. 4 
1,25 | Pobre i ay : 
} 
1.54 $234.8 Led4d.2 
t,25 $552.8 1584.5 { 
2.40 $553.5 


TABLE XIV. DEEP NOTCH FEA LINEAR RESULTS - NODAL 


os 4800 psi 


DISTANCE 
FROM NOTCH, in. Oo, spsi Oo. ,psi 


TABLE XV. DEEP NOTCH FEA LINEAR RESULTS - GAUSS 
ea 4800 psi 


i DISTANCE 


FROM NOTCH, in. Cy rpsi Oy psi 


MW. Se Sega. 4 eli. 
2 4,75 eas2,2 toe ae) 
7 

1.48 Base. 7 | Se atte 

125 a436.1 SSe4.7 

1.548 Si4t.3 Sadi. 8 

1.75 $520.1 w4f1.2 

2.80 4622.48 SESE. 
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TABLE XVI. SHALLOW NOTCH FEA NONLINEAR 60,000 lb LOAD ; 


DISTANCE NO LOAD RESIDUALS 


FROM NOTCH,in. oO. ,psi Oo. ,psi oy +psi oe *psi 


TABLE XVII. SHALLOW NOTCH FEA NONLINEAR 65,000 lb LOAD 


NO LOAD RESIDUALS 


DISTANCE 
FROM NOTCH,in. Fy *psi Oo. ,psi Ty *PSi S. *psi 


TABLE XVIII. SHALLOW NOTCH FEA NONLINEAR 70,0001b LOAD 


NO LOAD RESIDUALS 


DISTANCE 


FROM NOTCH,in. Tp psi o.. »PSi Oo, ,psi Oo. ,psi 


1.348 TS 941.5 2baeg.< t 2569, B38.8 
bits T4SR3P, 5 SHETALS Lao, 4 2o98.7 
2.48 SE1LT7. L34e8 6 1448.2 2207.8 
TABLE XIX. DEEP NOTCH FEA NONLINEAR 30,000 lb LOAD 


NO LOAD RESIDUALS 


DISTANCE 
FROM NOTCH, in. Sy »psi o. *PSi ao, ,psi o.. psi 


TABLE XX. DEEP NOTCH FEA NONLINEAR 35,000 lb LOAD 


NO LOAD RESIDUALS ' 


DISTANCE 
FROM NOTCH,in. oF +psi Oo. ,psi 0, »psi Oo. ,psi 


tee nbeeet 


| 
TABLE XXI. DEEP NOTCH FEA NONLINEAR 40,000 lb LOAD 
NO LOAD RESIDUALS 
iY DISTANCE 

‘ FROM NOTCH,in. J, spsi Oe psi Cy psi G2 PSL 
a,c 378.7 243335 a bec Pe | -2152.3 
i, 34 SeSre.€ Sesis. 4 ~e lia. 4 2377.4 

: rs e384. 8 248731 L468.2 143.1 
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1.28 S4n38, 5 $35.5 riQ.s { 

dea Se5it2.8 teerz.z 2a. 7 Hed. 3 

aw 21872, 2 15382,3 232.8 $78. 


89 


oO ioe SM 


» diately aby SON 


corp wee ss arama ene: 


Soe ne gate ais: 


TABLE XXII. RIGID - PERFECTLY - PLASTIC RESULTS 


ey = 73,000 psi 


DISTANCE 
FROM NOTCH ,in. oO, ,psi o. ,psi 


TABLE XXIII. EXPERIMENTAL DATA ) = 0.25 HOLE LINEAR LOADING 


o, = 10,749 psi 


DISTANCE 
FROM HOLE, in. 0, ,psi a +psi 
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TABLE XXIV. EXPERIMENTAL DATA SHALLOW NOTCH LINEAR LOADING 
| 
15,000 lb Load 

1 
DISTANCE 
FROM HOLE, in. ah Lik Oe +psi 
ema Tare) e 
| 12k Wat. 8 
‘ A, 258 2350 7 

ie iecn ase 

A. Sea Sa, SS 

W.625 {ara , 7 

TABLE XXV. EXPERIMENTAL DATA SHALLOW NOTCH 60,000 lb LOAD 
: NO LOAD RESIDUALS 
} DISTANCE 
; FROM NOTCH,in. 9, ,psi o epsi 0, psi o. «psi 
=~ : 
7 
1. ee Sride.d L3ers.4 
{ 
j 91 
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TABLE XXVI. EXPERIMENTAL DATA SHALLOW NOTCH 65,000 1b LOAD 


NO LOAD RESIDUALS 


DISTANCE 


FROM NOTCH,in. [% ‘Psi ry "psi S_ spsi o. psi 
= : 
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A, SPS B4863.2 | 1F435.1 


TABLE XXVII. EXPERIMENTAL DATA SHALLOW NOTCH 70,000 1b LOAD 


j NO LOAD RESIDUALS 
DISTANCE 
. FROM NOTCH ,in. Fy psi Oo. ,psi Ty psi o. ,psi 


4.258 S425, 3 {3713.2 “487 74.8 ~a7e.&@ 

| 
: 4.375 | a, a Lse1L5.8 ~ P4eer, “Lalla. él 
g.599 | 2524.6 f 16814.9 -EOHLT LA Petes 


W,boS Adda, 4 1r4di.n W247 950 -i5S 1.8 
4 


j Arse BH455.5 | Ss138.4 


B.375 | FESGS.F eae aa 


Se 
_ 
Dox] 
ams 
Dut) 
aap | 
ea 
Co 
ts 
L 
_ 
feo 
oa 
J 
™ 
Bon 
a 


ers EE Ly 


TABLE XXVIII. EXPERIMENTAL DATA DEEP NOTCH LINEAR LOADING 


15,000 1b Load 


DISTANCE 
FROM NOTCH,in. 


TABLE XXIX. EXPERIMENTAL DATA DEEP NOTCH 30,000 lb LOAD 


Elastic - Plastic 


DISTANCE 
FROM NOTCH,in. 
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APPENDIX B 


s LOCAL LEAST SQUARES SMOOTHING 


® corner Nodes 
o Midside Nodes 
+ 2x2 Gauss Points 


Two-Dimensional Isoparametric Element from ADINA [Ref. 4] 


The local smoothing expression from Hinton and Campbell 


[Ref. 19] in ADINA coordinates becomes 


oy Cc B B A Ty 
54 B A Cc B Ory 
~ = xX 
3 A B B Cc Srrr 
oy B Cc A B Sry 
_ V3 ahaa = 43 
where A= 1 + zs B=~- 5 and C=1- > 


With Tye Jos 33 and om representing the smoothed corner node 


7rr° Srrye and Try as the unsmoothed 


stresses at the Gauss integration points, this expression 


stresses and <7., 


can be written in an equivalent form. 


o A BO B o. 

7a A BOC : Seay 
ont . Cc B A B Sry 
om B Cc B A Cry 
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INCLASSIFIED 


The midside node stress values may be obtained by averaging 
the values at the associated corner nodes, since the distri- 
bution of the smoothed stresses is linear along the sides of 
the element. Smoothed stress values obtained by this least 
squares method should subsequently be averaged to obtain 


unique values at nodal points shared by adjacent elements. 
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